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Abstract 

We present the first attempt of a new method to compute the pion light-cone wave 
function (LCWF) on the lattice. We compute the matrix element between the pion 
and the vacuum of a non-local operator: the propagator of a "scalar quark" (named for 
short "squark" ) . A theoretical analysis shows that for some kinematical conditions (en- 
ergetic pion and hard squark) this matrix element depends dominantly on the LCWF 
$7t(u),m £ [0, 1]. On the lattice, the discretization of the parton momenta imposes 
further constraints on the pion momentum. The two-point Green functions made of 
squark-quark and squark-squark fields show hadron-like bound-state behaviour and 
verify the standard energy spectrum. We show some indications that during a short 
time, after being created, the system of the spectator quark and the squark behave 
like partons, before they form a hadron-like bound state. This short time is the place 
where the partonic wave function has to be looked for. 
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1 Introduction 



The light cone wave functions (LCWF) Q enter the calculation of a large variety 
of processes such as electroweak decays, diffractive processes, meson production in e + e~ 
and 77 annihilation, relativistic heavy ion collisions, heavy flavors and many others [Q]. 

The LCWF depends on a large momentum scale, n 2 , which is typically the momentum 
of the considered hadron P 2 in a physically well chosen reference frame (e.g. equal velocity 
frame for form factors, B rest frame for B decay, etc.). The pion wave function is expanded 
in terms of Fock states: 



k) = ai\qq) + a 2 \qqg) + a 3 \qqgg) H . (1) 

where the lowest Fock state \qq) describes the valence configuration which is dominant at 
large enough P 2 ||. Up to power corrections 0(Aq CD /P 2 ), the valence component \qq) 
is fully described by its leading twist amplitude. 

The leading twist amplitude has been proven to be describable in a very compact 
and frame independent way: the wave function & n (u) is defined by the following matrix 
element involving the 7r _ -meson and a light cone Wilson string 



(0\d(0)V 



exp(i f dr^) 

J X 



lnlbu{x)\-K (p)) x 2 =0 = -i Pfl f n / due %v *- x $„(u) . (2) 







The Wilson string in the square bracket ensures the gauge invariance of the l.h.s. of eq. 
@. The link between the first term in eq. (|l]) and 3V (u) will be dicussed in subsection 



2.1 



Let us notice here that eq. (Q) describes the LCWF a la Bethe Salpeter (BS), but, al- 
though the Bethe Salpeter framework differs significantly from the null-plane quantization 
approach, eq. (Q) exactly describes the dominant contribution to the pion wave function 
on a null plane. It is also useful to remind that the null-plane quantized wave function on 
a plane t + z = 0, is equal to the pion wave function quantized on t = constant, for a pion 
with a momentum P z = oo. In eq. (||) u denotes the longitudinal momentum fraction of 
the pion carried by the (valence) quark in the infinite momentum frame. The antiquark 
carries a fraction (1 — u). 

Let us insist, the pion wave function in QCD is an extremely complicated object, 
which cannot be reduced to the BS wave function on the light cone f\. However, in its 
infinite momentum frame, it simplifies dramatically in the following sense: the form factors 
depend only on the longitudinal wave function defined in eq. (|2|) while the transverse 
motion of quarks becomes irrelevant. For finite but large pion momenta the corrections 
are 0(Aq CD / P 2 ). Equivalently, for a quark and an antiquark lying almost on the same 
light line a corrective term 0(x 2 Aq CD ) has to be added to the l.h.s. of eq. (|2|) if this is 
not to be restricted to x 2 = 0. 



2 We use the expression "light cone" wave function according to a common habit although "null plane" 
wave function is more appropriate since the quantification surface is indeed a null plane. 
3 The light cone is a surface of zero measure in full space time. 
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Systematic expansions in inverse powers Aq CD /P 2 may be performed. But, even 
better, for each order in Aq CD /P 2 , perturbative QCD (pQCD) methods [|l],^],[5| allow the 
coefficients to be systematically expanded in powers of 1/ log(P^/Aq CD ). 

The dominant term in this perturbative expansion, i.e. the asymptotic form of the 
LCWF for very large fi 2 ~ P 2 reads: 

*?(u) = flu(l-u) (3) 

In this extreme limit, the shape of the wave function is totally given by pQCD, while 
the multiplicative constant f w in eq. @ contains all the relevant non-perturbative knowl- 
edge. The function (||) is corrected by terms which decrease only logarithmically when 
fi 2 — > oo. While the anomalous dimensions of these terms are computable from pQCD, 
their coefficients are only computable by non perturbative methods or to be taken from 
experiment. 

At lower jj?, when the 0(Aq CD /^ 2 ) power corrections can still be neglected but not 
the logarithmic 0(1/ log(/x 2 /AQ CD )) ones, the form of the wave function evolves away 
from eq The study of the LCWF in this range needs the use of non-perturbative 
methods. Most frequently one computes the LCWF via moments of the function ^^(u) 
as will be shortly described in the next paragraph. A well known example is the work by 
Chernyak and Zhitnitsky (CZ) || who used the QCD sum rules [] to determine the first 
two moments and obtained that at fj, = 1 GeV the shape of the pion wave function is 
completely different from its asymptotic form and it writes: 

$£ Z (u) = 120u(l - u)(u - 0.5) 2 (4) 

As can be seen from eqs. (|) and (|) there is a large difference between the two functions. 

Experimental measurements of the electromagnetic form factors of the pion were con- 
sidered to be the best way to study these wave functions |§. Recent model-dependent 
analyses of CLEO data on meson-photon transition form factors |), 1C] are consistent with 



the asymptotic wave function. A direct measurement [11] was carried out using data on 
diffractive dissociation of 500 GeV/c ir~ into di-jets from a platinum target at Fermilab 
experiment E791. The results show that the asymptotic wave function (||) describes the 
data well for fi 2 ~ 10 (GeV/c) 2 or more, although this interpretation is subject to some 



controversy [12 



On the theoretical side, a direct non perturbative measurement of the LCWF is badly 



wanted. There are only few attempts in that direction. The first method [13| is a lattice 
computation of moments of the LCWF 

M n = [ du u n 3> n {u) . (5) 
Jo 

which can be done computing the pion to vacuum matrix elements of local operators such 
as 

(TT-(p T )\d(0)^ 5 (iD^)...(iD^)u(0)\0) = -i/,KpK 1 -.# + "' 



4 These QCD sum rules for the first two moments of the pion twist-two distribution amplitude were 
recalculated in ref. Q resulting in a shape between the two extreme cases <!> as and $ cz . 
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where the dots at the end correspond to terms suppressed by powers of Aq CD /P^ (the 
same terms have been eliminated in the l.h.s. of Eq. (|2|) by means of the restriction 
x 2 = 0). The lattice discretization of the derivative operators in (|6|) is more and more 
tricky with higher moments, and their renormalization isn't easy either. 

It was therefore proposed in |H| to attempt a direct calculation of the LCWF from 
lattice QCD. One tries to "see" on the lattice the partonic constituents of the hadrons 
instead of the hadrons themselves. The idea is first to consider an energetic pion, which 
is supposed to have its partonic constituents "frozen" by Lorentz boost, and second to 
hit one of its quarks by giving it a large momentum in order to measure the perturbative 
part (small distance between the constituents) of the wave function. A scalar with the 
color content of quarks propagating from the hit quark to the spectator one insures gauge 
invar iance. 

In this paper we report the first and preliminary real attempt in that direction. In 
section |2| we explain the principle of the calculation and derive the basic formulae, with 
a particular care at establishing for which parameters of the run we may expect the 
subdominant contributions to the pion wave function to be under control. In section ^ 
we describe the lattice set-up which was used. In section ||| we present the results on the 
two-point Green functions. In section |H| we present the results on the three-point Green 
function and present the main analysis of our result. We believe that our results might 
provide some hint of a partonic behaviour. Finally, we discuss the relevance of our results 
in section |(| 

2 Principle of the calculation 

In this section we want to elaborate some theoretical tools necessary to prepare the 
direct lattice calculation of the LCWF. The issue is to reach some understanding of what 
to run on a lattice to measure the pion LCWF and to estimate the expected uncertainties. 
It is clearly impossible on a lattice to measure directly the matrix element in eq. (|2|) since 
obviously Euclidean metric has no light cone. The large momentum frame approach is 
more promising, with a standard continuation to imaginary time. We will then need to take 
into consideration the full pion wave function, assuming from QCD some general knowledge 
about it, and then consider under which conditions what is measured in the lattice depends 
dominantly on the LCWF, and if so, to estimate the subleading contributions. This will 
first be performed in Lorentz metric in an infinite volume. Later on we will take into 
account the Euclidean metric and the finite volume effects. 

2.1 Derivation of the basic formulae 

From now on, we will use the Light-cone gauge, where the path ordered operator 
V exp(z dr^A^) is equal to 1. Equation @ defines the pion Bethe & Salpeter wave 
function on the light cone, which has been extensively studied in literature since the 
pioneering work of Brodski and Lepage Q. It contains the leading contribution to the pion 
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wave function, the subleading pieces having been eliminated by the light cone condition 
x 2 = 0. We are aiming at a lattice investigation of this wave function. This will lead 
us (as already done in ref. ||14| ) to compute Fourier integrals of the wave function over 
the whole space and not only on the light cone. Therefore, the effect of subdominant 
contributions should be considered. Luckily, hadron properties, as derived from QCD 
asymptotic freedom, allow to control the approximation introduced when neglecting these 
subdominant contributions. 

Let us follow the standard Light-cone perturbation theory (LCPth) techniques pL 
We consider the first term in eq. (|l|), i.e. the valence u — d Fock state ^| for the 7r _ -meson 
wave function resulting from the quantification on the null-plane time i.e. x + = t + z = 

(v+n = Vo + (-)v z y. 

< 0\d(0) 7 + l5 u( x )|tt( p ) > x+=0 = -ip + U J due- iu ^ J |^ e^^ ud/7T (u, k ± ) 



= -ip + fn due~ luV ~^~ ip ud/n (u,x ± ) (6) 
J 

where the change of variable k + = up + has been performed, with < u < 1 since both the 
"+" components of quark (up + ) and antiquark ((1 — u)p + ) have to be positive (remember 
that components "+" of momenta have to be positive by definition) and where ipud/ir{ u i %±) 
is the partial Fourier transform (over kj_) of Vw/7r( u ! 

The previous matrix element depends on the light-cone three-momentum p = (p + ,p±) 
and its conjugated three- vector in configuration space, x = (x~,x±). For the sake of 
simplicity, we chose the frame where p± = 0, and hence p ■ x = p + x~ /2. The wave 
function ipud/Tr( u ^±) m Eq. represents the probability amplitude for finding two 
partons with momenta (up + ,k±) and (p + (l — u),—k±) respectively in the valence Fock 
state of the pion. This amplitude is normalized to 1, 



du 



o 



d 2 k\ 



j^Y^ud/A^ k -i) = J duip ud/7r (u,Q) = 1 , (7) 



as it immediately comes from requiring that < 0|(i7 + 75«|7r >= —ip + f n when the operator 
becomes local, i.e. when x = in Eq. (||). 

In eq. (|6|) we have only considered the 7 M component in the direction "+" of the pion 
momentum. The other directions j 1 - and 7~ can lead to matrix elements proportional 
respectively to p 1 - and p~ . In the pion rest frame all these components of the matrix 
elements should be of order Aq CD if we do not assume any restriction 0on x. This simply 

5 Strictly speaking, we retain only the dominant part of the valence Fock state, the one connected to 
vacuum via the axial current, the other contributions being suppressed. This suppression can be understood 
simply from the fact that the quarks in an energetic pion have dominantly the same helicity. 

6 Let us repeat that we are not allowed to restrict ourselves to small since we will perform Fourier 
transforms. 
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expresses that the size of the pion in its rest frame is O(Aqcd) in momentum space and 
0(1/Aqcd) in configuration space. Let us now consider a frame in which the pion has a 
very large p + . Then the matrix element considered in eq. (|^) is increased proportionally 
to the increase of p + , on the contrary x~ is decreased by the same ratio and the transverse 
components stay constant. For an "infinite momentum" pion we are left only with the 
contribution proportional to the pion momentum. This is a first indication that in our 
analysis we will have to concentrate on energetic pions. 

Equation. (||) is a definition of the wave function i^ud/-n{ u i k±). It only depends on the 
quantities u (the fraction of pion's momentum carried longitudinally by one parton) and 
k±; it is frame- independent for longitudinal boosts. In order to establish the connection 
with eq. @, we now put x 2 = (i.e. x± = provided that we quantized on the light-cone 
time x + = 0) in eq. (|6|). If we take x± = in eq. @ and compare the result with eq. (g) 
we see that [] 



/d 2 k± - ~ 

j2^2^ad/n(u,k±) = ^ud/n(u,0) . (8) 

To clarify the physical picture let us now compare in the Light-cone gauge the l.h.s of 
eq. @ unrestricted to x 2 = (the full BS equation) and the l.h.s of eq. @. They only 
differ by the null plan constraint x + = 0. This constraint is generated by requiring that 
the pion carries a large momentum. Indeed p~ = m 2 /(p z + E n ) appears to be powerly 
suppressed. This suppression of p~ implies that + p~x + ~ (unless x + is 

unnaturally large). If one assumes the absence of sudden changes when x + moves away 
from 0, one may replace p + x- by p x in eq. @ which now reads: 

<0|J(0)7 M 75«(x)|7T-(p) >= -ip^U I due- iup - x ip ild/n (u,x ± ) , (9) 

J o 

If we add the physical input that the wave function extends typically to transverse mo- 
menta of the order of Aqcd , we get from < u < 1 the picture that the valence constituents 
of the pion move essentially in the same direction than the pion itself at a velocity close 
to 1. In other words, due to asymptotic freedom, the constituents do not like to have a 
very large virtuality and the only way for almost massless quarks to build up the energy 
and momentum of the almost massless pion is to move in the same direction, i.e. to have 

Eg + E q ~ \k q \ + \kq\ ~ E n ~ \k q + k q \ 



From now on we shall follow the method in [14] and we will replace the gauge invari- 
ance restoring operator V exp(i dr^A^), by another one, which is easier to continue 
analytically to euclidean time: the scalar coloured propagator, 

S{0; x) = ' ' 



T> 2 — m 2 s + ie —d 2 — m| + ie 

d 4 k -ifc-fO-o;) « ( 1Q ) 



(2^) 4 k 2 - ml + ie 



Remember that the exponential in brackets in eq. (H) is equal to 1 in our gauge. 
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where ms is a mass parameter, assumed to be small or zero to mimic a massless parton. 
In eq. (|io|), when replacing T> 2 by d 2 we have bluntly neglected the coupling to gluons. 
This has been done in order to simplify the argument which will follow and is justified if 
we assume the scalar object to be "hard" and hence to behave mainly as a parton. Still, 
a careful study of the effect of radiative corrections is strongly needed. This replacement 
looses the gauge invariance of the 1/T> 2 operator. This is difficult to avoid: if the light 
cone wave function (||) is gauge invariant, the more general ones, eq. (|9|) are not. The loss 
of gauge invariance is here the price we pay to present the argument which will follow. 
Needless to say, the real lattice calculations have been performed in a gauge invariant way. 

We can thus write 



e-""* <0|J(Ofr„76S(0;aOii(i)|*(p) > 
This is supposed to be valid for all x so that we can integrate over x and obtain: 



% \ d A x e~ iq x < 0\d(0)^ l5 S(0;x)u(x)\7r(p) > 



o 



1 1 2vr (2vr) 2 kj, - k 2 - m 2 + te ' 



l r a2 



i(uEK-\J (up z +q z ) 2 + (q ± +k ± ) 2 +m^ t 



d k± e 

= -P/ifn I du I — . VW/ttO, (12) 

Jo J [Z7T) 2^(up z + q z ) 2 + (q ± + k ± ) 2 +m 2 

where = 0, xo = —t (t < 0), kn = (ko,k z ) and again p± = 0. The r.h.s. of the latter 
line derives from integrating the former's over kn and changing variables (q — k)± — > —k±. 

At this stage let us return to the physical understanding of the wave function ipud/ir( u i ~k±) 
already briefly considered above. The quarks have a small probability of being far off shell 
and ipud/Tr( u ^ k±) vanishes when k\ becomes large 0. In practice, k 2 1 _'tp a( j j /-K{ u ^i-) ^ as 
k\ —> oo 0. Therefore, this suppression for large k\ allows one to expand in powers of 
the transverse components, provided that E$ S> Aqcd; Aqcd being a natural hadronic 
energy scale bounding the transverse momentum carried by the partons and 



E s = J(up z + q z ) 2 + q{+m 2 s . (13) 



s Perturbative analysis indicates that hadronic wave functions do not decrease quickly enough as k\ — » 
oo to avoid the appearance of infinities. The Pion qq-wave function falls off roughly as [JjJ, and the 

resulting UV logarithmic divergence is the origin of the scale dependence of the wave function. For the 
sake of simplicity this point shall be deliberately overlooked in our formal derivation. 
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We then get: 



/ / d'.v f-"'- e <0\d(0)^ l5 S(0;x)u(x)\ir(p) > 



/•l e i(uE n -E s )t 



+ 



d 2 k ± 



2E% 



1 \ lPud/Tr(u, ~k±) 



+ ... , 



(14) 



It is easy to see that the second term inside the bracket (last line) is formally O (Aqcd /Eg), 
provided that t <C Es/Aq CT) and t <C Eg /(Aqcd\q±\)- This second term is negligible as 
long as, and this is the general situation, Eg ~ [\up z + q z \ 2 + Aqcd 0- However, 

when \up z + g z | ~ Aqcd for some values of u and when <7±^Aqcd ; i-e. when p and (f are 
back to back the expansion in eq. (14) breaks down as Eg is not larger than Aqcd 
any longer. In another language, giving a large transverse kick to the pion generates a 
hard gluon exchange between quarks which selects the perturbative component of the pion 
wave function, the so-called "small pion", which is what we want to measure. Indeed in 



the FNAL experiment E791 [11], the LCWF is observed via jets which have rather large 
transverse momenta. Let us now summarize. 



Conditions for a partonic signal: [CI] In order to reach some knowledge on the 
light-cone wave function ^(u) from the lattice calculation of the l.h.s of eq. (|i~4|), the 
following conditions are required beyond the general large pion momentum constraint i.e. 
Pz » Aqcd : t <C £s/A QCD , t <C Eg / (Aq r>\q±\) and Eg > A QC d for all u. This 
generally implies cos m i n <cos6 pq for some cos m j n significantly greater than -1. 

2.2 Consequences of discrete partonic momenta. 

Let us now consider a finite parallelepipedic volume with periodic boundary conditions 
(torus). As is well known, the momenta components can only take the form 

= -r n » ( 15 ) 

where are integers and is the length in the direction [i. This is obviously valid also 
for partonic momenta FJ. Thus in the formulae of subsection 2.1 all integrals over J n X du 



have to be replaced by discrete sums over the values of u such that up^ verifies eq. 

Here comes immediately a problem. Let us assume for one moment that the com- 
ponents of p^ are all or 2ir/L^. Then only the values u = 0, 1 are allowed. In any 

9 remember that ms is small. 

10 Strictly speaking p and q could be back to back as long as \q z \ — \p z \ 3> Aqcd. 
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For other values the amplitudes are canceled by destructive inteferences. 
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model the LCWF which is proportional to u(l — it), (^) and (H), vanish for these values. 
The expected dominant behaviour at large momentum is vanishing in this case, and only 
subdominant effects can be observed. 

The simplest situation, the only one considered from now on, is when the pion mo- 
menta are aligned along one of the lattice spatial directions /x. To allow values of u that 
scan the domain of variation [0, 1] densely R enough to provide a fair description of the 
LCWF we should have : 

2vr , , 

P„ > -T- ( 16 ) 

This condition [C2] has to be added to the set of conditions [CI] summarized at the 
end of subsection |2.1| in the case of infinite volume. Clearly this new one is not equivalent 
to the former ones since this one does depend on and disappears smoothly in the large 
volume limit. 

2.3 Strategy for lattice calculations 

Following the method of 14] we compute on the lattice the three point Green function 



F»(p,q;t)= / d 3 yd 3 xe-^e- l vy < 0|P 5 (y, ^(x, t)S(x,t; 0) W d(0) |0 > e^^-*).(17) 



When all the conditions [CI] and [C2] are satisfied, and after performing a Wick rota- 
tion to euclidean metric, the l.h.s. of equation (|l7|) verifies approximately the following 
proportionality in terms of the LCWF: 

F»(p, <?; t) cx Pfi U Y, $7r W (18) 

Hi 

where the ^ extends over all values < Ui < 1 such that Uip^ * L/(2ir) are integers, p 
is the momentum of the pion generated by the interpolating field Ps(y) = d{y)^§u{y), q is 



a momentum given to one valence quark of the pion. E${u) is defined in (13). We have 
assumed < t < t n . The g^C^-*) takes into account the propagation of the pion between 
t and t n . Of course t n — t has been assumed to be large enough to eliminate the excited 
pseudoscalar states. 

Eq. (^) may be understood in a simple way : the time evolution between and t 
is the product of the propagators of two "partons" , one scalar parton of energy E$ with 
a propagator proportional to e~ Est /(2Es) and the spectator quark of energy (1 — u)E v . 
The scalar parton has the color quantum numbers of a quark. For convenience let us call 
it a squark although it has obviously nothing to do with supersymmetry. The three-point 
Green function in eq. ( |l7|) could also be used to estimate the form factor for the transition 

12 The dominant contribution to the LCWF is only possible when all the n M ,/x = 1,3 are or have a 
common divisor, and at least one n M is larger than 1. 
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between a pion and squark-quark bound state (which we call a pionino, tt, to follow on 
the same metaphoric nomenclature). In such a case we would take t large enough for the 
ground state pionino to dominate: 

F^(p,q;t) « e~ Eit (19) 

t— »oo 

where is the pionino energy. For small t, on the contrary, the excited states should 
add up coherently in a complicated manner. The analysis presented in the preceding 
subsection seems to indicate that this should boil down to a rather simple partonic-like 
picture. In other words we expect a kind of hadron-parton duality to be at work for small 
t which should allow a partonic reading of our data. At this stage it is clear that we need 
to study, beyond the three-point function in the l.h.s of eq. (|l~7|), the two point function 
corresponding to the pionino interpolating field. 

An additional comment concerns the squark mass. In the preceding formulae we have 
written a squark mass as a free parameter. In order to gain the richest possible 
information on the pion wave function, the renormalized squark mass has to be as light as 
possible. How to perform this ? We have chosen an approach based on an analogy with 
QCD hadrons. We will vary the bare squark mass down to when the algorithm to compute 
the squark propagator stops converging, which we take as an indication of possible zero 
modes. 

Finally, all things considered, we will have to make a systematic study of the spectrum 
of all the colorless bound states constituted by quarks and squarks. It will turn out that 
in the quenched approximation nice exponential behaviours do indeed appear, signaling 
the existence of pioninos and squark-squark bound states (see Fig. [l]), and furthermore, 
for non vanishing momenta, they follow the relativistic spectral law E = y m 2 + p 2 (see 
Fig. |2j) , or if one prefers the lattice one (see eq. ^2|) below), which is not distinguishable 
from the former within our statistical errors. 

2.4 Symmetries 

Before turning to the actual calculation, it is useful to summarize which among the 
two- and three-point Green functions we intend to compute should vanish because of 
QCD's discrete symmetries. 

In general all Green functions we consider are real in configuration space. Therefore 
they are real in momentum space if parity-even and imaginary if parity-odd. See table |j. 

3 Lattice set-up 

We consider a 16 3 x 40 lattice at /3 = 6.0 in the quenched approximation. The quarks 
are computed with the clover action with the coefficient c sw = 1.769. We have used for 
the spectator quark two values of the bare mass parameters: k = 0.1333 and 0.1339, and 
for the active one k = 0.1339. 
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operator 


7 matrices 


parity 


real/im 


time reversal 


vanishes at p = 


squark-squark 


t 


+ 


real 


+ 


no 


squark-quark 


t 


+ 


real 


+ 


no 


squark-quark 


7o 


+ 


real 




no 


squark-quark 


7i 




imag 


+ 


yes 


quark-quark 


75 - 75 


+ 


real 


+ 


no 


quark-quark 


7o75 - 75 


+ 


real 




no 


quark-quark 


7*75 - 75 




imag 


+ 


yes 


three-point 


7o75 - 75 


+ 


real 




no 


three-point 


7i75 - 75 




imag 


+ 


yes 



Table 1: This table shows the symmetry properties of the Green functions. By three-point we 
mean the Green function F^lp, q) defined in eq. (p"7|). The second column refers to the 7-matrices 
in the Green function. For squark-quark, only one 7-matrix is traced with the quark propagator. 
In the other cases we indicate the matrices on both ends of the quark propagators. The third 
column refers to the spatial parity of the Green function. The time reversal refers to the symmetry 
when t — y — t (and — > — t w in the the three point case). We thus learn, for example, that the 
three point with 7075 — 75 vanishes at t = if t w = t max /2. 



The squark propagator D(x,0) verifies the equation 



,y-fi + ~~ A)$ 



x,y+p, 



D(y,0)=5 xfi (20) 



We compute the squark propagator with the bare mass parameter k$ = 0.1428, 0.1430, 
0.1431. Above k$ = 0.1431 the convergence of the inverter becomes very long, which we 
take as a sign that we are close to the massless squark. 

In each case we have run 100 configurations. The errors are computed according to 
the jackknife method. The pion interpolating field P5 is inserted at t n = 16. This has been 
chosen so that the direct signal at small t is not significantly perturbed by the signal which 
has looped around via the end of the lattice : 40-16 has eight time-intervals more than 16. 
This is an important precaution. Indeed from table [l] we learn that the three-point Green 
function with 7075 inserted at t = is odd for time reversal. If was taken in the middle 
of the lattice, t = 20, it would have resulted a vanishing of this three-point Green function 
for t = 0. Since we are interested in small values of t such a vanishing of the signal would 
have made the analysis impossible. 

For the study of the two and three point Green functions we have run with the following 
values for the pion three momentum: 

^ = (0,0,0); (1,0,0); (1,1,0); (2, 0,0) (21) 

In practice, however, the vanishing momentum does not produce a pion describable by a 
light-cone wave function. The momentum (1,0,0) (1,1,0) will not be useful since in these 
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pionino two point Green function 



pionino two point Green ftmctio 
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10. 20 30. 40. 



Figure 1: Two point Green functions in logarithmic plots for pionino and squark-squark states. 
As an example we present the lightest states, i.e. n — 0.1339 and ks = 0.1431. 

cases only the values u = 0, 1 are allowed by the discretization and the LCWF vanishes 
for these values. However, they are kept in the analysis for a comparison of the results 
obtained from (1,0,0) (1,1,0) with the ones from (2,0,0), which might be interpreted as 
partonic signal. 

Concerning we have run a large number of momenta, with components ranging 
from — to ^ but again too large momenta are too noisy. We will detail later the 
momentum configurations which are considered in the analysis. 

As already explained, we hope to catch the partonic signal at small t. In practice we 
have concentrated on the region t = 0, 4 as we will see later. It leaves — t > 12 which 
should be enough to isolate the pion and it leaves some space to look for plateaus. 



4 Two point Green functions 

We have shown in fig. |l| six examples of new two point Green functions for momenta 
p = (0, 0, 0), p = 27r/L(l, 0, 0) and p = 2ir/L(2, 0, 0) respectively. It is seen that these two 
point Green functions do indeed behave as if the quark-squark and squark-squark states 
were hadron-like bound state. 
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momentum 





1 


1.4 


2 


pion qiqi -75,75 


0.42(2) 


0.62(3) 


0.70(4) 


0.61(13) 


pion q x q x - 7o75,7o75 


0.41(2) 


0.60(2) 


0.69(3) 


0.89(7) 


pion q 2 qi -75, 75 


0.38(2) 


0.60(3) 


0.66(4) 


0.35(16) 


pion q 2 q2 - 75,75 


0.34(2) 


0.58(4) 


0.61(5) 


0.09(22) 


pion q 2 q 2 - 7o75,7o75 


0.34(2) 


0.56(3) 


0.62(4) 


0.84(10) 


rho q x q x -7^7* 


0.62(1) 


0.76(2) 


0.96(3) 


1.02(6) 


rho q 2 q 2 - 7i>7« 


0.60(2) 


0.71(3) 


0.98(5) 


0.95(9) 


pionino q\Sl — 70 


0.59(1) 


0.71(1) 


0.81(1) 


0.98(3) 


pionino q x Sl — 11 


0.55(1) 


0.67(1) 


0.77(2) 


0.91(5) 


pionino gi52 — 70 


0.54(1) 


0.67(1) 


0.77(2) 


0.95(3) 


pionino q±S2 — 11 


0.51(1) 


0.63(2) 


0.74(2) 


0.88(6) 


pionino qi 53 — 70 


0.51(1) 


0.65(2) 


0.76(2) 


0.93(3) 


pionino q\S3 — 11 


0.48(2) 


0.60(2) 


0.72(2) 


0.86(7) 


pionino q 2 51 — 70 


0.57(1) 


0.70(1) 


0.79(1) 


0.98(3) 


pionino 52 51 — 1 


0.53(1) 


0.64(2) 


0.74(2) 


0.91(7) 


pionino q 2 S2 — 70 


0.52(1) 


0.66(2) 


0.76(2) 


0.94(3) 


pionino q 2 S2 — 11 


0.48(1) 


0.60(2) 


0.71(3) 


0.87(9) 


pionino q 2 S3 — 70 


0.49(2) 


0.63(2) 


0.74(2) 


0.92(4) 


pionino q 2 S3 — i 


0.45(2) 


0.57(2) 


0.69(3) 


0.85(10) 


squark-squark 5151 


0.59(2) 


0.70(2) 


0.80(2) 


0.93(5) 


squark-squark 5252 


0.50(2) 


0.61(2) 


0.74(3) 


0.83(7) 


squark-squark 5353 


0.44(2) 


0.56(3) 


0.72(4) 


0.74(8) 



Table 2: Energies of the various bound states in units of a -1 (for (3 = 6.0, a -1 ~ 2.0 GeV). 
The symbols (71,(72 represent respectively n = 0.1333,0.1339 for quarks; 5*1,52, 53 respectively 
ks — 0.1428,0.1430,0.1431 for scalars. The momentum norms are given in units of 2tt/L. We 
indicate the 7 matrices used in the meson interpolating fields. 

We present the results for the energies of the bound states in table ^. In fig [2] we present 
some checks of the spectral law E = \J m 2 + p 2 . The latticized free boson dispersion 
relation 

sinh 2 (£/2) = sinh 2 (m/2) + ^ sin 2 (p M /2)) (22) 

does not significantly differ from the continuum one within our errors. For momentum 
4tt/L the quark-quark states are in some cases meaningless due to the noise. It is surprising 
that the non conventional states present a better signal for this large momentum. 

Of course the main lesson of this analysis is that the non-conventional bound states, 
pioninos and squark-squark do behave exactly as real hadrons. We are not in a position 
to discuss the theoretical implications of this fact neither make any statement about the 
existence of such bound states in a non-super symmetric extension of QCD. 
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Light meson energy 



Light squark-squark state 




0.0 .20 .40 .60 

P 



Light pionino state 



Heavy pionino state 



e=0.1339 
K.=0.1431 



K=0.1333 
K„=0.1428 



- i 



- i 
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P 



Figure 2: Energy of the bound states as a function of the momentum in lattice units. The dots 
correspond to the continuum formula E — y/ m 2 + p 2 , the mass being taken as the central value of 
the zero momentum energy. Three first plots are with k = 0.1339, ks — 0.1431, the last one with 
k = 0.1333, Kg — 0.1428. 



The lowest bare squark mass considered is k,$ = 0.1431. When ks is varied slightly 
above 0.1431, the scalar inverter does not converge anymore. This squark is coded 53 in 
table ^ and we see that the corresponding squark-squark bound state rest mass is about 
0.44 in lattice unit, i.e. about 900 MeV (a -1 ~ 2 GeV for j3 = 6.0), not far from the rho 
meson mass. It is rewarding that the mass of this squark-squark bound state is rather 
light, as if the squark with an approximately vanishing renormalized mass did indeed 
produce rather light bound states 0. We feel encouraged to treat indeed this squark as a 
light parton as will be done soon. 



We do not know of any symmetry which would impose a pion-like massless state for massless squarks. 
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5 Three-point functions 



With our set of momenta, only the momentum p n = 2tt/L (2,0,0) gives a non- 
vanishing [j &-k(u) for discrete u = 1/2. Thus we will focus our analysis on the latter 
momentum although we have studied the full set of momenta p n , with a set of momenta 
q to be discussed later. We have only considered the time component F°(p,q;t). 

Our analysis of the data follows from section 2^ . To test whether eq. or eq. ([l9| ) 



has some relevance for our data we will consider whether the following quantities: 



F°(p,q;t) 

and 



p -(E 7r /2+E s )t 
P fn ^ MV2) 



(23) 



fVwit)^'] 1 (24) 



are constant in time for some time interval. 

Before that, it is instructive to have a look at the numerators F°(p,q;t). In figure ^ 
we have plotted as an illustrative example the three-point function for p n = 2ir/L (2, 0, 0) 
and various vectors q. We observe a very striking feature akin to an oscillating behavior. 
We do not claim to understand fully this shape. However since, in section a rationale 
was elaborated to describe the expected partonic behavior which may show up at small 
time, we will focus from now on on this time interval. 

The very rapid drop observed at small time, i.e. t E [0, 3 — 4] is present for all values of 
q. We will test the hypothesis that this rapid drop is due to a partonic signal assuming 
that the hadronic behaviour sets in for larger times. The typical shape in fig. ||] might 
suggest a negative interference between the small time regime and the later one, leading 
to a vanishing amplitude around t = 4. We do not understand the origin of the latter, 
which is beyond the scope of this work focused on the small-time drop. It is noticeable 
that the statistical errors for this time range are small enough to exhibit a signal while 
the two-point function for the corresponding pion propagation time and the same pion 
momentum is extremely noisy. 



Searching for plateaus at small times: A plateau of eq. ( p3|) would indicate a 
partonic-like behavior, while a plateau of eq. ( |2"4| ) would sign a pionino. We will compute 
E v and E% from the measured pion and pionino rest masses, see tab. ||, and the formula 
E = y m 2 + p 2 . We prefer this to the direct use of the measured energies for non zero 
momentum, reported in tab. ^ because the latter are noisier than the rest masses for 
p= f (2,0,0). 

The energy Es has been taken via eq. (|l3l ) assuming two possible masses ms for ks = 
0.1431. As already mentioned, for n$ > 0.1431 the calculation of D(x, 0) from eq. ( |20| ) 

14 Notice that the CZ wave function (^) vanishes for u — 1/2, its study needs even larger momenta and 
will not be discussed in our analysis. 
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fails indicating the presence of small eigenvalues, i.e. that mg is small. Besides considering 
a massless scalar parton (ms = 0), we have also considered the value ms = 0.22 in lattice 
units, which corresponds to the scalar-scalar bound state mass (divided by two). It would 
be tempting to fit mg from the results yielding the flattest plateau, but it turned out to 
be too difficult to disentangle the effect of ms on the plateau from other effects which will 
be discussed later. 

In fig. |] we show two examples of ratios corresponding to eqs. (p3|) (left) and ( |24| ) 
(right) at small time. In the light of the discussion in subsection |2.1| , we have chosen as il- 
lustrative the following kinematics: Lq/(2n) = (—2, —1, —1) and Lq/(2ir) = (—1, —1, —2), 
both for Lp/(2ir) = (2,0,0). It is clearly seen that that the plots to the right (|24j) are 
utterly incompatible with a plateau, thus discarding a pionino interpretation at small 
time. The plots to the left might show some indication of plateaus but they deserve some 
discussion. The signal decreases from a maximum at time to reach a value compatible 
with at a time 3 — 4. This happens not only for these two examples but is a general 
pattern for all the kinematics considered. This cancellation has already been seen on the 
numerators of eq. (^3|) in figs. ||. We have argued that it is motivated by destructive 
interferences which generate an overdecreasing of the numerators in eq. (^) with respect 
to the denominators. The signals vanish as soon as t = 3 — 4, restricting the range where 
plateaus might be seen to a very short time interval ^] around t = 0. 

Anyhow, the most restrictive of the constraints relative to t, summarized at the end 



of subsection 2T, i.e. t <C E$ /(Aqcd\q±\), amounts, for a massless scalar parton, with 



our lattice set-up and the value u = 1/2, to the condition 



a 



-1 



(Px/2 + qx? + q\ 



t<£- -* : ; 5 (25) 

A QCD I'LL | 

where for Aqcd we have taken a typical quark transverse momentum of 400 MeV within 
a hadron. This constraint does not allow to use larger time domains than the one just 
discussed. 

We will now go on confronting the slopes on this small time interval to the theoretical 
prediction of a plateau for eq. (p3|), postponing the maybe more convincing comparative 
study of the values of F°(p, q; 0). 

We perform a systematic study over a larger set of three point Green functions defined 
such that: Lp/2n = (2,0,0), (L<f/2vr) 2 < 4 and (L (q+p)/2ir) 2 < 6. These limitations on 
the norm of the momenta are meant to avoid too noisy results. On the other hand, the 
constraint Eg 3> Aqcd (see subsection |2.l| ) translates into the lower bound: 



L \(Vx , _ \ 2 , 2 



, , +fe +gi>i (26) 

7T V \ Z ' 



One may worry about contact terms or other lattice artifatcs that might spoil the analysis around 
t = 0, this will be discussed in the conclusion. 
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For this set of data we measure the slope of the ratios in eqs. ( p3| ) and (24) for the 
time intervals t = 0, 3 and t = 0, 4. For the latter range, the results are presented in fig. 
^ the ratios of eqs. ( p3| ) and (p4|) are presented for commodity as a function of the cosine 
of the angle between p and q, which we will from now on refer to as cosO pq . 

We have eliminated from the analysis the data with Lqj '(2tt) = (—1,0,0) for which 
the scalar parton is at rest (p x /2 + q x = 0) and thus violates the condition eq. (|26|) . 
The data with white circles on the plots correspond to Lq/(2n) = (—1,0,-1) which is 
marginal for both conditions eqs. (25,^6). It could be noticed that the back-to-back points 
Lqj '(2tt) = (—2,0,0) do not raise problems as a result of the discretization of partonic 
momenta. Indeed, since u = 1/2, up + q never vanishes contrarily to the continuum case 
discussed in section 2T. More generally, the majority of the points with cos6 pq close to -1 
are not excluded for the same reason. 

Comparing both plots in fig. [|, it is evident that the partonic slopes (left) are much 
closer to zero than the hadronic ones (right). Nevertheless, the partonic slopes show a 
general tendency to be negative (see tab. ||) which can be traced back to the vanishing 
around t = 3 — 4. The white circle points show a lesser improvement of the partonic data 
as compared to the hadronic ones as conjectured just above. 



model 


time slice 


xVcLo.f 


average slope 


pionino 


0-4 


4.1 


-0.56(18) 


partons ms = 


0-4 


1.9 


-0.26(13) 


partons ms = 0.22 


0-4 


0.92 


-0.23(13) 


pionino 


0-3 


8.3 


-0.82(9) 


partons ms = 


0-3 


3.7 


-0.39(13) 


partons ms = 0.22 


0-3 


2.2 


-0.36(13) 



Table 3: Average slopes (and % 2 /d.o.f for a vanishing slope) of the expression appearing in eqs. 

) and ( pi[ ) for two time slices and two parton masses. It is seen that the parton mass does not 
play a very important role. The difference between the two time slices is due to the zero of F° 
discussed in the text. 



The slopes given in table ^ are the averages over our set of momenta q (excluding the 
momentum corresponding to the white circle). We have kept the mass of the scalar parton 
between and half the mass of the scalar-scalar bound state (see tab. |2|). The resulting 
slopes do not depend significantly on the latter mass. It can also be seen that the slopes 
are quite similar for time-slices [0,3] and [0,4]. 



Comparing three-point functions at t 

features of the partonic behaviour: 
i) the exponential time evolution 



= 0: Equation ( |l8|) predicts two main 
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ii) the following amplitude at t = 



The begining of this section was devoted to the time evolution. Let us now focus on the 
amplitude eq. (p7j). 

The plot in fig. |6| shows for our set of momenta q the product Es(l/2)F° (p, q\ 0) which 
is expected to be constant from eq. (|?]). Eg is computed from eq. ( |T^ ) with a massless 
scalar parton. The plotted ratio is indeed strikingly constant: the y^/d.o.f. for the fit 
to a constant ratio is 0.22. This expected constancy of a large set of numbers, which are 
significantly different from zero, yields an amazing support to a partonic interpretation 
of these data. We cannot figure out any other explanation for this feature. Indeed, 
one might fear that the observed constancy of Es(l/2)F° (p, q; 0) is simply due to some 
contact term producing a q independent F°(p,q;0) combined with a small dependence 
of Eg{l/2) on q. To consider this we have tried a fit with F°(p,q;0) = constant, which 
gives x 2 /d-o.f. = 0.72 , larger than the previously found 0.22, although still smaller than 
1. We would thus rather believe, in agreement with the partonic interpretation, that the 
small variation of F°(p,q;0) is a consequence of the constancy of Es(l/2)F°(p, q; 0) and 
a small variation of Es(l/2). As a check, we have tested the constancy of F°(p,q;0) 
for p = 27r/L(l, 0, 0), which is not expect to follow eq. ( |27| ) while contact terms have 
no reason to be absent 0. We find x 2 /d.o.f. = 2.7 which further supports the partonic 
interpretation of the constancy Es(l/2)F°(p, q; 0) for p = 2ir/L(2, 0, 0). 



6 Discussion and Conclusion 

We have performed the first tentative application of a new proposal [14] to compute 
the pion LCWF. This proposal was to compute the pion to vacuum matrix element of a 
non-local operator, namely the propagator of a scalar particle which has the color quantum 
numbers of a quark. For convenience, we call it a "squark" . This resulting matrix element 
is gauge invariant. To exhibit the partonic structure of the pion a large momentum q is 
added to the scalar propagator. 

We have shown that, provided the pion has a large enough momentum p, provided 
that the squark has a large enough energy, and provided the propagation time of the 



scalar object is short enough (end of subsection 2.1), the above mentioned matrix element 
is dominated by a contribution from the pion LCWF. A measure of this matrix element 
can then provide informations on the LCWF. 

A necessary first step is the computation of the two point Green functions of quark 
(squark) - quark (squark) bound states. The new states, which contain at least one 
squark, show a behaviour quite similar to standard hadrons, they show nice exponential 



16 We did not check the constancy of Es(l/2)F°(p, q; 0) in this case since u = 1/2 is forbidden in the 
case p = 2tt/L(1,0,0). 
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time dependence, fig. |l], they verify Einstein spectral law, fig. 0, and the masses decrease 
with increasing k$ i.e. decreasing squark bare mass. 

We have then analyzed the three point Green functions for a large set of pion momenta 
p and transfers q. The scalar parton has a momentum up + q where u G [0, 1] is the fraction 
of pion momentum carried by the active quark. The discretization due to the finite volume 
implies a discretization of the fraction u. In our set, only the momentum p = ^(2,0,0) 
allows for (where the LCWF vanishes), namely u = 1/2. 

We focused the analysis on small times (t E [0, 3 — 4]) according to the formulae fl23| ) 
and ( p4[) which express respectively the hypothesis of a partonic behaviour of the squark 
and the spectator quark during this small time interval, or, on the contrary, the hypothesis 
of a precocious confinement of the squark and the spectator quark into a hadronic-like 
bound state. The correct hypothesis should show up as a plateau in time. 

Our data clearly favor the partonic behaviour at small time: the observed rapid drop 
of the Green function is expected from a partonic picture, while a hadronic picture predicts 
slower decrease. The analysis is however made delicate due to an observed vanishing of 
the Green function around t = 3 — 4 which might be due to a destructive interference. The 
resulting analysis domain is very short and close to zero. This might induce the objection 
that we cannot disentangle our signals from lattice artifacts such as contact terms, etc. 

Nevertheless, a second series of tests has confirmed our feeling that a real partonic 
signal shows up: all the Green functions at t = for our set of values of q verify the 
prediction, eq. (p7|), of the partonic model (up to one unknown constant) in an amazing 
manner. It is difficult to figure out how a lattice artifact could mimic this behaviour for 
so many data. 

This work aimed mainly at testing the viability of this program. We believe that the 
answer is positive. The fact that we could argue rather firmly that we see a partonic signal 
is encouraging, obtained on a small lattice, with a rather large lattice spacing, and "large" 
momenta which are indeed not so large ! 

In order to progress we first need to settle the question of possible lattice artifacts. 
To that aim, it would be necessary to change the lattice parameters and mainly a and 
to run a larger set of momenta. This would furthermore allow to reach other values of u 
than 1/2 and provide an idea about the shape of the LCWF. This program implies the use 
of a larger volume which would also hopefully reduce the noise of large momenta Green 
function. 

A recent work by S. Dalley based on a Hamiltonian formulation of QCD on a lattice fll5|] 
presents an interesting analysis of the LCWF. This new method is very promising although 
it presents some difficulties as stated by the author. It is of course too early to perform 
a detailed comparison of the Lagrangian formulation used here and the Hamiltonian one. 
Both need to be followed. 
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Figure 3: We plot F°(p,q;t), normalized by a constant (divided by the pion propagator with 
p = (An/ L, 0,0) from the fixed time t v to 0) versus the running time for momenta indicated on 
the plots using the lightest quarks (re = 0.1339) and the lightest "squark" (res = 0.1431).. 
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Figure 4: Ratios of eqs. ( p3| ) (left) and ( |24| ) (right) for momenta indicated on the plots using the 
lightest quarks (k = 0.1339) and the lightest "squark" (ks = 0.1431). 
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Figure 5: Slope of the ratios on the time interval t=0,4 for formulae 
for different values of q and for p = ^(200) with a massless scalar parton 
the cosine of the angle between vectors p and q. 
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Figure 6: Values of Es(l/2)F°{p, q; 0) normalized as fig. || for p — (An/L, 0, 0) and our full set of 
q (labeled from 1 to 12 on the horizontal axis). The data show the expected constancy around the 
average represented by the horizontal line. 
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